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ABSTRACT 

Peculiar velocities of galaxies hosting Type la supernovae generate a significant sys- 
tematic effect in deriving the dark energy equation of state w, at level of a few percent. 
Here we illustrate how the peculiar velocity effect in SN la data can be turned from 
a "systematic" into a probe of cosmological parameters. We assume a flat A-Cold 
Dark Matter model {w — —1) and use low and high redshift SN la data to derive 
simultaneously three distinct estimates of the matter density which appear in the 
problem: from the geometry, from the dynamics and from the shape of the matter 
power spectrum. We find that each of the three f^m's agree with the canonical value 
Q.m = 0.25 to within Icr, for reasonably assumed fluctuation amplitude and Hubble 
parameter. This is consistent with the standard cosmological scenario for both the 
geometry and the growth of structure. For fixed — 0.25 for all three l^m's, we 
constrain 7 — 0.72 ± 0.2f in the growth factor ^rn{z)'^ , so we cannot currently distin- 
guish between standard Einstein gravity and predictions from some modified gravity 
models. Future surveys of thousands of SN la, or inclusion of peculiar velocity data, 
could significantly improve the above tests. 
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1 INTRODUCTION 

The observed present acceleration of the universe was first 
confirmed a decade ago by two separate groups using Type 
la su pernovae (SN la, iPerlmutter et al.l Il999l : [r iess et al.l 
Il998h . SN la are one of a number of probes needed to obtain 
tighter constraints on dark energy equation of state, includ- 
ing any possible time evolution. This will require surveys of 
thousands of supernovae out to high redshifts to accurately 
measure their luminosity distances from which parameters 
describing the dark energy can be inferred. To achieve the 
desired constraints on dark energy, in particular a few per- 
cent constraint on the equation of state parameter w, the 
supernovae will have to be accurately calibrated. It is there- 
fore vital that this calibration is done accurately, and it is 
the low redsh ift supernovae which are vital to achieve this, 
for details see lAldering etall (|2002l ). At low redshift the su- 
pernovae distances have little or no dependence on the cos- 
mological parameters such as Q,m, and the dark energy 
equation of state w. They do however put a tight constraint 
on a combination of what is essentially the calibrated mag- 
nitude zeropoint (A/) and the Hubble constant Ho, whereas 
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for the high redshift supernovae there is a strong degen- 
eracy between M, Hq and the cosmological parameters of 
interest. Figure [T] illustrates the importance of the low red- 
shift supernovae in a nchoring the Hubbl e diagram. It shows 
the gold sample from lRiess et al.l (|2007l ) constraints on Qrn. 
and w with and without supernovae with redshifts less than 
0.1. One can see that without the low redshift supernovae 
(blue/light contours, using 146 SN la) the constraints blow 
up significantly compared to the full gold sample (red/dark 
contours, using 182 SN la). There are several sources of sys- 
tematic error which affect the calibration of the zeropoint, 
for example dust extinction, l uminosity evoluti on, weak lens- 
ing, and Malmquist bias (see lKim et al.l[2004l . for more de- 
tails). This type of error is not decreased by having a large 
number of supernovae and will necessarily come to dominate 
the error budget. The systematic errors mentioned above 
have long been discussed in the literature and are not con- 
sidered in this Letter. There is a source of error which is 
unique in the fact that it affects only the low redshift "cal- 
ibrating" supernovae, their peculiar motions relative to the 
Hubble flow. 

Previous authors have set about using SN la to quan- 
tify the degradation of dark energy errors due to pecu- 
liar motions or use them to trace the peculiar velocity 
field itself in a vaxiety of ways. Three distinct approaches 
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to this have been discu s sed re cently in the hterature. In 
iNeiU. Hudson, fc ConlevI ([2007^ different flow models based 
on the IRAS PSCz survey ([iBraiichini ct al. 1999|) were used 
to "correct" the luminosity distances by the known peculiar 
velocities before fitting them for the cosmological parame- 
ters of interest. They find the potential systematic error in 
w caused by ignoring peculiar velocities is of the order of 4 
per cent, i.e. quite significant. 

iRadburn-Smith. Lucev. fc Hudson! (|2004l ) compared pe- 
culiar velocities from 98 local supernovae with the gravity 
field predicted from IRAS. In HauEb0lle ot al. (2007) an an- 
gular expansion of the radial velocity field was used to probe 
the local dipole and quadrupole of the velocity field at three 
different distances. They found that the dipole is consistent 
with galaxy surveys (e.g. lErdogdu et al.l [20061 ) at the same 
Hubble flow depths. 

The third and s o mewh a t different method is ut - 
lised by iHui fc Greene! (|2006| ): ICoorav fc Caldwelll l|2006! ): 
[Gordon. Land fc" Slosaj (|2007l ) who take a covariance matrix 
approach. From the fluctuation in the luminosity distance 
induced by the peculiar motions , (see Hui fc Greene! 12006! : 
iPvne fc Birkinshaw! |2004 1 1996! : ISasaki! 1 1987! . for deriva- 
tions), a covariance matrix for the resulting errors in the lu- 
minosity distance (or similarly the apparent magnitude) can 
be calculated. The covariance matrix depends on cosmolog- 
ical parameters which describe the growth and distribution 
of structure. In addition to the peculiar velocity effect this 
is due to gravitational lensing effect, which is important for 
reds hifts larger than 1, and we shall ignore it in this Letter. 

ICoorav fc Caldwell! (|2006l ) found that peculiar velocities 
of the low redshift supernovae may prevent measurement of 
w to better than 10 percent, and diminish the resolution 
of the time de rivative o f w pr ojected for planned surveys. 
[Gordon. Land fc Slosa^ (|2007l ) used the covariance matrix 
approach on current data, showing the changing constraints 
on (Tg, fim and w depending on the exact redshift range 
of the SN la sample and whether the full covariance was 
included or not. They also apply the analysis to forecasting 
constraints for future surveys. 

Here we unify the analysis of SN la data to study si- 
multaneously fits for the expansion of the universe and the 
growth of structure. There is plenty of discussion on the 
possibility that the accelerated expansion of the universe is 
caus ed by a modification of general relativity on large s cales 
(e.g. iDurrer fc Maartendliooj : iHuterer fc Linder!l2007l . and 
references therein). By measuring the growth of structure, 
which directly effects the observed peculiar velocity field, 
information is gained to differentiate between the two sce- 
narios. 

The rest of the Letter is organised as follows. In Section 
[2] we describe the SN la sample used in this Letter. In Sec- 
tion 13] we describe the theory underlying SN la analysis in 
cosmology and the effect of the velocity field. 



2 DATA 

We analyse both nea r by su pernovae [z < 0.12) from 
|jha. Riess. fc Kirshnerl ([2003) and high redshift supernovae 
(z < 0.176) from a sample c ompiled by iDavis et al.l 



C^^gold sample 
V gold sample with z>l 




Figure 1. Assuming a flat universe with a constant equation 
of state w, these are 1 and 2 a likelihood contours showing the 
constraints in the Om-io plane for the Riess et al 2007 supernova 
data. The red/dark contours are when using the full gold sample 
(182 SN la), the blue/light contours are when using the gold 
sample but removing all supernovae with redshifts less than 0.1 
(146 SN la) 




Figure 2. The 1 and 2 a contours on S'Ud the parameter 

7 from using all 271 low and high redshift SNe. Values of other 
parameters include = 1 and h = 0.7, o-„ = 300kms~^, am = 
0.1 and the other two f2m's=0.25. 



l|2007r and 

IWood-Vasev et al.! l|2007! ). iDavis et all l|2007h combined 



I2OO7!) which includes data from iRiess et al 



the data from the two samples by normalising to the 
low redshift supernova e they had in common. Following 
|jha. Riess. fc Kirshneij ([2003) 9 supernovae are excluded 
from the low redshift set, those that are unsuitable due 
to bad lightcurve fits. This includes supernovae with their 
first observation more than 20 days after maximum light, 
those that are hosted in galaxies with excessive extinction 
{Ay > 2.0 mag) and one outlier (SN1999e), which appears 
to have an extremel y large peculiar velocity. This l eaves 124 
supernovae from the ljha. Riess. fc Kirshneil (|2007l ) data set 
in the redshift range z £ [0.0023,0.12], and median red- 
shift z = 0.017. The overlapping SN la in the two data 
sets were used to est imate a smal l normalising offset to the 
magnitudes from the lDavis et al.! ((2007. ) data set (the extra 
magnitude error is negligibly small). The same procedure 
was used bv lDavis et al.! (|2007l ) in normalising the two high 
redshift data sets. After eliminating duplicated SN la, our 
combined data set has 271 SNe with z G [0.0023, 1.76], and 
z = 0.29. 
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3 METHODOLOGY 

We describe here how we utihse the SN la dataset described 
in Section [2] to estimate cosmological parameters by includ- 
ing the peculiar velocity covariance. 



3.1 Covariance matrix approach 

The luminosity distance is defined as 



(1) 



where JT is the observed flux of the supernova and C is its 
intrinsic luminosity. The apparent magnitude m of a super- 
nova at redshift z depends on the luminosity distance as 
follows 



m{z) ^blog^o Dl{z) - 51ogio(//o) + Af + 25 



(2) 



where M is the magnitude zeropoint, and Dl is defined 
without the Hubble constant as Dl = Hodi in kms~^. The 
equation above ignores the additional terms which involve 
applying dust corrections, K corrections etc. For a flat uni- 
verse containing a matter component and a dark energy 
component with a constant equation of state, the luminosity 
distance can be written as 



DLiz) 



c{l + z) 



dz' 
eW) 



(3) 



E{z') = + + + (4) 

If the universe was truly homogeneous and isotropic (FRW) 
this would be the end of the story, the observed Dl would be 
described accurately by Eq. |31 However peculiar velocities 
have the effect of perturbing the luminosity distance 



5Dl 
Dl 



Vr 
C 



c(l 



H{z)dL{z) 



(5) 



where Vr is the radial peculiar velocit y of the supern o va an d 
H(z) is the Hubble paramet er . See lHui fc Gr eene (200(d); 
iBonvm. Durrer. fc Gasoarinil (|2006l ): IPvne fc B irkin shaw. 



l|2004 h Sasakil (1983) for a derivation. We emphasize that 
in the right hand side of the above equation di is for an 
unperturbed FRW universe, derived at a perturbed redshift 
z. 

Therefore the covariance of the perturbation in Dl, 
5Dl/Dl for a pair i,j is given by 



and 



1 - 



c{l + zf 
H{z)dL{z) 



1 - 



cjl + zY 
H{z)dL{z) 



(6) 



(vriVrj) ~ = cos9iCos9i'^^^{r) + sinOi sinOj'^ ±{r) (7) 

is the linear theory r adial peculiar velocity correlation func - 
tion, (jCorskil ITosj : iGroth. Juszkiewicz. fc Ostrikej Il989^ . 
The angles in Eq.[7]are defined by cos9x = fx - f and the di- 
agonal elements are given by Eq.[9]below. The and 
^x(r) can be calculated from the matter power spectrum 
using linear theory 



*ll,x(r) = Z)'(20^'fe) 



27r2 



Bi,j_{kr)dk 



(8) 




Figure 3. Constraints on Ojjf a^nd flm" from using all 271 low 
and high redshift SNe. In both panels 0^°™ has been fixed at 
0.25, o-„ at 300kms~^ and Cm at 0.1. In the top panel rig has been 
fixed at 1.0, and the red (dashed) contours, black (dark) contours 
and blue (light) contours are the results when h = 0.6, h = 0.7 
and h = 0.8 respectively. In the bottom panel h has been fixed 
at 0.7, and the red (dashed) contours, black (dark) contours and 
blue (light) contours are the results when rig = 0.9, rig = I and 
Ug = 1.1 respectively. The black (dark) contours in both panels 



are exactly the same. The grey line indicates where 



■■ n 



dyn 



where P{k) is the matter power spectrum, Bi_ = jo{x)/x 
and £?! | — jg (x) and jg , jg are the first and second derivative 
of the zeroth order spherical Bessel functions respectively 
and D'{z) is the derivative of the growth function at redshift 
z. D'{z) is a function of the Hubble parameter and Q.m- See 
Section 13.31 for further discussion. The auto correlation is 
given by 



i--\D'\z.) j^dk. 



(9) 



We can therefore calculate Cij for a pair of supernovae at Zi 
and Zj respectively given a set of cosmological parameters. 

In the above equations the growth factor is calcu- 
lated exactly numerically. More insight to the dependence 
on Q.m is given by the commonly used approximation for 

the growt h factor J_ = d\n5/d\na « Q.rn{zy , where 

7 « 0.6 (|Peeblesl ll98C ), with little dep endence on the 
cosmological const ant (Lahav et al.l 199ll). and a slight 
dependence on w (^W ang fc Steinhardd 199^ ). Recent re- 
fined calculations predic t 7 = 0.55 for the concordance 
model, and 7 = 0.69 l|Linder fc CahnI |2007| ) for a par- 
ticular modified gravity model, P GP braneworld gravity 
l|Dvali. Gabadadze. fc Porratil I2OO0I ) . though this is just an 
example of many possible modified gravity models. Below 
we shall constrain 7 from the SN la data. 
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3.2 Likelihood analysis 

To find the set of cosmological parameters ©max = [^i---6'n] 
that best fit the data we find the set that maximise the 
likehhood function. Assuming that the data and the obser- 
vational errors are Gaussian random fields the likelihood 
function can be written as 




(10) 



where Di is defined as A = {Of" ~ Dl{z)) /Dl{z) and S 
is the cov ariance matrix including the o bservational noise. 
Following iGordon. Land fc Slosad (|2007l ) we write this as 



(11) 



where ai is the standard uncorrelated error given by 
'ln(lO) 



(o-m + (A*r''')^)+ 



^(1 + ^)' I l±(i2) 

H{z)dL(z) / . C2 ^ ' 



where is often set to 300kms^^ and is included to ac- 
count for nonlinear contributions to gi, (w hich is derived 
only in linear theorv. [Silberman et ahlboOll ). and the veloc- 
ity of the SN within the host galaxy. Here Om is the intrinsic 
magnitude scatter and ^'''^^ is the error from the light curve 
fitting. 



3.3 The Three Faces of r2„ 



From the equations in Section [3T1] it can easily be seen that 
C/j is a function of ilm through 

(i) the geometry of the universe, from H(z) 
and also dhiz) in Eq.O For a flat universe fim"™ — 1 — Ha, 
with no dependence on other cosmological parameters, 
j-jgeom xaost strongly constrained however by the high 
redshift SNe through Eqs.[3]and|l 

(ii) fJjlf"; the growth of structure in the universe, 
from D'{z) in Eqs. [S] and O or the equivalent growth 
factor paramaterization / = [Q'l^"'{z)Y' .We note a strong 
degeneracy through the product asQm{z)'' where 



rimiz ^ 0){1 + zf 



(13) 



n„,{z = o)(i + z)3 + i- nm{z = 0) 

for a flat universe, our results can be scaled accordingly. 

(iii) Clm'- the matter power spectrum P{k) in Eqs. [8] 
and [9] It is wel known that the shape of the power spectrum 
depend on the product F — Q,^h, with some degeneracy 
with e.g. the spectral index Us , as and baryon and neutrino 
mass densities Qb and Q^- 

Please note that Eq.'s|8] and [9] contain all of the low redshift 
Qrn terms. If the ACDM model of the universe is correct 
then when varying each of these "faces" of fim separately 
the results should be consistent with each other. If not this 
suggests that the ACDM model is inconsistent and the data 
may favour a model which changes the theory of general 
relativity on large scales or other dark energy models. 



Table 1. Results for each "face" of Q^n under different param- 
eter combinations. Because we do not marginalise the errors are 
small, but it is still useful to look at the relative errors for the three 
Cm's. Columns A, C and D set the other two dm = 0.25, whereas 
column B marginalises over them. In column A and B the "nui- 
sance" parameters are ay = SOOkms"'^ and am = 0.1. Column 
C assumes the "nuisance" parameters are small {a^ = 200kms~^ 
and am = 0.08). Column D assumes the "nuisance" parameters 
are large (ay = 400kms~^ and am = 0.12). Each column has set 
the power spectrum parameters [h Us eg] = [0.7 1.0 0.8]. 





A 


B 

0.25lHi 


C 


D 

24+'''°i 


r-fdyn 


'-'•^0-0.14 










43+0-46 


0.38l°;« 


on+0.43 

u.oo_g 23 


0.48l°:g 



4 RESULTS 

In the following analysis we assume a flat ACDM universe 
with a dark energy equation of state w = —1. To gain an 
insight for the effect of varying the other cosmological pa- 
rameters, namely Ho, ris, and the "nuisance" parameters dm 
and ay we do not marginalise over them but present the re- 
sults at some choice values for these parameters. The effect 
of marginalising over cr„ and am degrades the error on °™ 
by about 10 percent, the error on Q,^ changes negligibly and 
the error on Qm" degrades by about 30 percent. The larger 
degradation of the error on Qm" is because of its strong 
degeneracy with ay. Other parameters are fixed as follows; 
Qb ~ 0.04, = and as ~ 0.8. For clarity most contours 
have only the la confidence level, and where appropriate the 
assumed values of other parameters are stated. 

Table [1] shows a set of results for each "face" of Qm 
under different parameter combinations. This table shows 
that the different Qm's are consistent with the canonical 
value of 0.25 to within 1-a. One can see the degeneracy 
direction of the errors {ay, am) and Qm", discussed below. 
See table caption for an explanation of the columns. 

Figure [2] shows the constraints on Qm" and the param- 
eter 7, described at the end of Section 13.11 The contours 
shown are the 1 and 2 a. For Qf^"^ = 0.25 the 68 percent 
confidence constraint on 7 is 7 = 0.72 ± 0.21 consistent 
with both Einstein gravity (7 — 0.55) and DGP gravity 
(7 = 0.69). 

Figure |3] shows the effect of h and Us on the Q,m-Qm" 
contours. As expected there is only a strong infiuence on 
Qm- One can also see the degeneracy direction of h and 
Us in the power spectrum, shown by the indistinguishable 
differences between the red (dashed) contours in both pan- 
els. The positions of the red (dashed) contours show that 
decreasing h by roughly 10 percent is equivalent to increas- 
ing Us also by roughly 10 percent. This is also shown by 
the blue (light) contours. The contours for Q,^-Q^°"^ and 
Qgeom_^dyn showu here since the same combinations 

of h and as plotted in Figure[3]have a nearly negligible ef- 
fect on the contour positions. Under different permutations 
of h — 0.6,0.7,0.8 and = 0.9, 1.0, 1.1 the contours only 
shift in either the f2JJf or Qm" direction and even then there 
is only a maximum shift of Aflm of 0.05. 
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Figure 4. Constraints on all the Qrn pairs from using all 271 low and high redshift SNe. From the left the first panel shows Q^-Q^", 
the second panel f2m-f2m°™, third panel f2^^°'"-C^" and the fourth panel O^^-cr^. The black (dark), red (dashed) and blue (light) 
contours in the first 3 panels correspond to (t„ equal to 250kms~^, 300kms~^ and 350kms~^ respectively. Each time ns = 1, h = 0.7, 
= 0.1, erg = 0.8 and in the first 3 panels the third Clm is kept fixed at 0.25 and just the l-a contour is shown. In the fourth panel 
ns = 1, h = 0.7, (J,n = 0.1, Qm = ^m°"^ = 0.25 and both the 1 and 2-a contours are shown. The grey line indicates where = Qm- 



Finally Figure |3] shows the effect of (t„ on all the contour 
pairs. It has the largest effect on iljjf" and very little effect 
on $1^^°'". This is because the diagonal elements of Cfj (see 
Eq. O are approximately proportional to (fij^")'^''' and to 
cr^. Therefore the best-fit OSS'" decreases as a„ increases. 



5 CONCLUSIONS 

We have presented in this paper a unified approach for prob- 
ing both the expansion of the universe and the growth of 
structure with SN la data and to test the consistency of the 
ACDM model. We utilised the SNIa data to derive three 
distinct estimates of the matter density Qm which appear 
in the problem: from the geometry, from the dynamics and 
from the shape of the matter power spectrum. We found 
that each of them agrees with canonical value flm = 0.25 
to within la. We note we are restricting our discussion to 
ACDM, if we allow w to vary our constraints on will 
weaken. We also constrained 7 in the growth factor r2m(z)^ 
and found for il^ = 0.25, 7 = 0.72 ± 0.21. This value of 
7 is consistent with both concordance and some proposed 
modified gravity models. 

Current and future SN la surveys such as SN Factory, 
GAIA and Skymapper (for low redshift), SDSS-II (for in- 
termediate redshift) and DES, Pan-STARRS, LSST, DUNE 
and SNAP (for high redshi ft) will generate s a mples of thou- 
sands of SN la (e.g. Alb recht etahl I2OO6I : iPeacock et"aLl 
I2OO6I . for overviews) . Large samples of low redshift SNe will 
greatly improve our constraints on and and the 

high redshift SNe on ilf^""^. Utilising galaxy peculiar veloc- 
ity data (using D„ — a and TuUy- Fisher distance indicators) 
will also provide improvement on the fiJJ^' and flfi^" con- 
straints. Our approach can also be generalised for a range 
of other cosmological parameters and exotic models of dark 
energy and gravity. 



6 ACKNOWLEDGMENTS 

We are grateful to Kate Land for providing us with infor- 
mation about SN la data sets, and to Sarah Bridle, Josh 
Frieman, Chris Gordon, Saurabh Jha, John Marriner and 
Jochen Weller for useful conversations. AA acknowledges 
the receipt of a STFC studentship and OL acknowledges a 
Royal Society Wolfson Research Merit Award. Both AA and 



OL thank Fermilab and the Kavli Institute for Cosmologi- 
cal Physics at Chicago University for their hospitality. We 
thank the anonymous referee for their detailed comments. 



REFERENCES 

Albrecht A., et al., 2006, astro, |arXiv:astro- ph/0609591| 
Aldering C, et al., 2002, SPIE, 4836, 61 
Bonvin C, Durrer R., Gasparini M. A., 2006, PhRvD, 73, 
023523 

Branchini E., et al., 1999, MNRAS, 308, 1 
Bridle S. L., Crittenden R., Melchiorri A., Hobson M. P., 
Kneissl R., Lasenby A. N., 2002, MNRAS, 335, 1193 
Cooray A., CaldweU R. R., 2006, PhRvD, 73, 103002 
Davis T. M., et al., 2007, ApJ, 666, 716 
Durrer R., Maartens R., 2008, GReCr, 40, 301 
Dvah G., Gabadadze G., Porrati M., 2000, PhLB, 485, 208 
Erdogdu P., et al., 2006, MNRAS, 368, 1515 
Gordon C, Land K., Slosar A., 2007, PhRvL, 99, 081301 
Gorski K., 1988, ApJ, 332, L7 

Groth E. J., Juszkiewicz R., Ostriker J. P., 1989, ApJ, 346, 
558 

E[augb0lle T., Hannestad S., Thomsen B., Fynbo J., SoUer- 

man J., Jha S., 2007, ApJ, 661, 650 
Hui L., Greene P. B., 2006, PhRvD, 73, 123526 
Huterer D., Linder E. V., 2007, PhRvD, 75, 023519 
Kim A. G., Linder E. V., Miquel R., Mostek N., 2004, MN- 
RAS, 347, 909 

Jelinsky P., SNAP Collaboration, 2006, AAS, 209, #98.09 
Jha S., Riess A. G., Kirshner R. P., 2007, ApJ, 659, 122 
Kolb E. W., Matarrese S., Notari A., Riotto A., 2005, 

PhRvD, 71, 023524 
Lahav O., Lilje P. B., Primack J. R., Rees M. J., 1991, 

MNRAS, 251, 128 
Linder E. V., Cahn R. N., 2007, APh, 28, 481 
NeiU J. D., Hudson M. J., Conley A., 2007, ApJ, 661, L123 
Peacock J. A., Schneider P., Efstathiou G., Ellis J. R., 

Leibundgut B., Lilly S. J., Mellier Y., 2006, astro, 

arXiv:astro-ph/061090 
Peebles P. J. E., 1980, The Large Scale Structure of the 

Universe, Princeton University Press 
Peel A., Knox L., 2003, NuPhS, 124, 83 
Perlmutter S., et al., 1999, ApJ, 517, 565 
Pyne T., Birkinshaw M., 1996, ApJ, 458, 46 
Pyne T., Birkinshaw M., 2004, MNRAS, 348, 581 



© 0000 RAS, MNRAS 000, 000-000 



6 Abate and Lahav 



Radburn-Smith D. J., Lucey J. R., Hudson M. J., 2004, 

MNRAS, 355, 1378 
Riess A. G., et al., 1998, AJ, 116, 1009 
Ricss A. G., ct al, 2007, ApJ, 659, 98 
Sasaki M., 1987, MNRAS, 228, 653 

Silberman L., Dekel A., Eldar A., Zehavi I., 2001, ApJ, 557, 
102 

Sugiura N., Sugiyama N., Sasaki M., 1999, PThPh, 101, 
903 

Wang L., Steinhardt P. J., 1998, ApJ, 508, 483 
Wood-Vasey W. M., et al., 2002, AAS, 34, 1205 
Wood-Vasey W. M., et al., 2007, ApJ, 666, 694 



© 0000 RAS, MNRAS 000, 000-000 



